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1. Introduction 
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Abstract. We use the recent supernova data set from the ESSENCE collaboration 
combined with data from the Supernova Legacy Survey and nearby supernovae to test 
the DGP brane world model and its generalisations. Combination of this data with 
I a flatness prior and the position of the peak of the CMB disfavours the DGP model 

ly-^ ' slightly. Inclusion of the baryon acoustic peak from the Sloan Digital Sky Survey 

0^ ■ increase the tension of the DGP model with the data, although it is not clear how 

self consistent this procedure would be without a re-analysis of the survey data in the 

■ framework of the DGP cosmology. Generalisations of the DGP model are tested and 
f — , constraints on relevant parameters obtained. 

o 

^ l' E-mail: sararydQphysto . se, malcOcern. ch, arielOphysto . se 

o 

Most of the energy in the universe seems to be dark, and many different observations 
^ . suggest that only a minority of this energy can be dark matter, the rest apparently 
being some kind of energy component with negative pressure. 

The simplest explanation for this dark energy component is a cosmological constant 
and people have attempted to explain the smallness of this energy density and hence 
its relatively recent dominance using anthropic arguments [1] which may concur with 
modern predictions from string theory [2]. Other explanations exist however, including 
a class of theories known as quintessence where the observed acceleration of the universe 
results from the stress energy of some rolling scalar field which only comes to dominance 
recently [31 H] . 

A different approach is to ask whether General Relativity breaks down at large 
distances, in other words to modify the left-hand/curvature side of the Einstein 
equations rather than the right-hand/stress-energy side. One of the most studied 
examples of such a scenario is the DGP brane world model [5], [6], [7]. This paper aims to 
compare the predictions of this theory and its generalisations with the latest cosmological 
data. 
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In the DGP model, gravity is trapped on a four dimensional brane world at short 
distances, but is able to propagate into a higher dimensional space at large distances. 
The Lagrangian for the model is {h = c = 1, Mpi = (SvrG)^^) 

Due to the different mass scales M and Mpi, gravity propagates differently on the brane 
and in the bulk. The effect of gravitational leakage into the bulk will only appear at 
large distances {r > L = M'^J2M^). 

The tt Friedman equation in this theory takes the form 

2 H pm . , 

where e = ±l. Itise = +1 that gives the late-time accelerating solutions that could 
explain dark energy. The theory also predicts very small but potentially detectable 
modifications to the earth- moon distance [8]. There are also changes in structure 
formation which could give constraints upon the theory [9] (see also |T0]). 

Recently the self- accelerating branch of this theory has come under theoretical 
attack because it seems to possess ghost-like instabilities [11] while one of the original 
authors has argued that the calculational regime within which the instabilities are found 
is not valid . In this work, we will put this issue to one side and proceed to compare 
the model's predictions with the latest data. Even if it turned out that the DGP model 
was theoretically suspect, it is interesting to see how different expansion histories may 
or may not be ruled out by the developing suite of data. 

To test theories, we would like to make detailed comparisons with astronomical 
observations. There are basically two ways of doing this. As already mentioned, one is 
to look at the growth of perturbations in these models and see how they compare with 
perturbations in ACDM and in the observed galaxy correlation function [9]. Another 
way is to simply look at the space-time geometry of the universe on the largest length 
and time scales available in order to reconstruct the expansion history of the universe 
and compare it to the solutions of the modified gravity models. This latter approach is 
the subject of this paper. 

Studies regarding the comparison between data and the DGP model based upon 
expansion history are to date contradictory. The first paper on the subject was presented 
by two of the current authors [13] using the 2005 data release from the SuperNova Legacy 
Survey (SNLS) [H]. In that work it was argued that the SNLS data combined with the 
position of the baryon acoustic peak in the Sloan Digital Sky Survey were less consistent 
with a fiat universe in the DGP model than in the ACDM model. These conclusions 
were backed up when Maartens and Majoretto performed the same test using the SNLS 
data, the baryon acoustic peak and the CMB shift parameter |15j . 

Later, a new "gold" data set of supernovae was released p^. Analysis of this data 
and the CMB shift parameter suggested that a fiat universe was completely consistent 
with the DGP model [ITJ, a conclusion the present authors can confirm, also upon 
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addition of the data from the baryon acoustic peak. It is interesting to note that it has 
been suggested that there may be an inconsistency between the different parts of this 
gold data set [T8 ]. 

The two sets of data come from different instruments - the SNLS supernovae are 
detected using a combination of imaging at the Canada-France-Hawaii Telescope and 
spectroscopic studies on large ground based instruments, in particular Gemini, VLT 
and Keck. The Riess et al. 07 gold sample contains supernovae from the Supernova 
Cosmology Project, SNLS, the High-Z Team, the GOODS transient survey, and includes 
21 new supernovae at extremely high redshift obtained with the Hubble Space Telescope. 
In all cases, low redshift supernovae from independent surveys have been included in 
the data sets. 

Furthermore the supernova magnitudes are obtained from the two different sets of 
supernova data using two different methods, the SALT algorithm which is used by the 
SNLS group fits supernovae based upon their light curve and their colour as does the 
MLCS method favoured by Riess et al. , the most recent version being called MLCS2k2 
|20] . The parameters in these two algorithms are obtained by training on low redshift 
supernovae before they are used to obtain the magnitudes of high redshift supernovae. 
SALT and MLCS2k2 differ both in the details of the brightness-shape relation correction 
and in how extinction/colour corrections are applied. 

Recently the ESSENCE group has released a set of 60 supernovae at intermediate 
to high redshifts [21]. This data has been combined with the SNLS data and data 
at low redshift to form a new data set spanning a large redshift range in detail. The 
supernovae in this data set have been analysed using both methods - SALT and MLCS. 
In this paper, we combine this new data with the CMB and baryon oscillation data in 
order to see if the DGP model is favoured or disfavoured. 

While the model as it stands has one extra space dimension, it is possible to imagine 
generalisations of the model with higher dimensional bulks, although because of the 
curvature singularity an infinitely thin 4D brane would create in more than 5 space- 
time dimensions, the theory would need to be regularised at the location of the brane 
[22] . We will present a parametrisation of these higher dimensional generalisations as 
has been done in previous work [23l [24l US] and we will introduce an extra degree of 
freedom which may go some way to modeling the regularisations. 

2. Observational status of Expansion History. 

Observations of high-redshift Type la supernovae have been used for over a decade to 
map the expansion history of the universe [25l[Ml[23[28l|29l[30l[3ll|32l[33l|3l[35l 
ESI [an |3Hl [391 [HI iOl Uni Elfl. These high-z surveys, combined with low redshift SN 
sample to anchor the Hubble diagram [l2l [l3l [H], provided the first direct evidence 

I In fact, the very first detection of a high-z Type la SN, SN1988U at z=0.31, was done almost twenty 
years ago by [41j . However, the data on this SN was scarce. Thus, this object is normaUy not included 
in the compilation of SNela distances. 
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that the universe expands at an accelerated rate. Cross-cutting techniques involving 
the anisotropies in the cosmic microwave radiation (CMB) [l5l |36[ HTl HH HHj, mass 
density estimates from X-ray observations of clusters [50] and more recently weak lensing 
[51], baryon acoustic oscillations (BAO) [521 [53] and the integrated Sachs- Wolfe effect 
[HU [55] are apparently continuing to confirm our current understanding of the universe, 
namely it being spatially flat, as predicted by inflation, and the current expansion being 
dominated by dark energy. 

The ESSENCE group has analysed supernovae at intermediate redshifts (0.1 < 
z < 0.8) using both the MLCS2k2 and the SALT methods to obtain the supernova 
magnitudes from the light curves and spectra. When using the MLCS2k2, the group 
adopt what they refer to as the 'glosz' prior to model Ay - the extinction in the V 
band of the supernova in the host galaxy [21j. They have also fitted higher redshift 
supernovae (0.11 < z < 1.1) from the SNLS survey [14] and a set of nearby supernovae 
{z < 0.11) again using both the SALT and MLCS2k2 method to obtain magnitudes of 
supernovae over a large redshift range. 

Measurements of the brightness of Type la supernovae as a function of redshift are 
sensitive to the cosmological model via the integration over expansion history in the 
expression for the luminosity distance (c = 1) 

,, = ^i^s(vmr4^) (3) 



ffoVM y Jo -^w 

where the function S{x) is defined as sin(a;) for fi^ < {closed Universe), sinh(a;) for 
> {open Universe) and S{x) = x, and the factor is removed for the flat 

Universe. The parameter E{z) = H{z)/Hq. 

In our analysis we fit the cosmological model to the supernova data using the 
luminosity distance above but we also fit to the baryon acoustic peak detected in the 
SDSS Luminous Red Galaxy survey (LRG) of Eisenstein el al (2005) [52l HQ], which 
constrains the following combination of parameters 



M 



E{zi) 



1 ^ ( r— dz 

-.S x/Wk 



^ivW \ 'io E{z 



0.472 ± 0.017, (4) 



where Zi = 0.35 and S and Qk are defined as in Eq.([3]). The quoted uncertainty 
corresponds to one standard deviation, where a Gaussian probability distribution has 
been assumed. 

There is some debate [T71[58] as to whether or not one should use the baryon acoustic 
peak to constrain models of dark energy which behave differently to a cosmological 
constant, for two reasons. The first is that the reconstruction from redshift space to 
co-moving space required to accurately identify the position of the acoustic peak has 
been done assuming a constant equation of state [52]. While one would expect the 
change in the position of the acoustic peak in an alternative dark energy model where 
the equation of state is a function of redshift to be small, the correction is difficult to 
quantify without detailed study for each model in question. Secondly, in modified models 
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of gravity, one would expect structure formation to proceed differently [H]. Although 
at the first approximation this would not change the physical co-moving size of the 
acoustic peak feature in the correlation function, such an effect might create systematic 
distortions in its reconstruction from redshift space. 

For these reasons we will constrain the models with and without the SDSS baryon 
acoustic peak, we leave it to the reader to decide if they want to pay attention to the 
BAO constraints on parameter space or not. 

Finally, a cleaner measure of geometry is the CMB shift parameter [561 US] (see 
however [57]) - the expansion history of the universe has to be such that the observed 
position of the CMB peak corresponds to the physical horizon size at last scattering. The 
photons used to measure this angular size have passed through the integrated geometry 
created by the particular dark energy model in question. The angular size of the first 
peak of the CMB as measured by the WMAP 3 year data constrains the ratio [591 EO] 



^0 ^y^) 

which can then be applied as another cut to parameter space for each model. 

Having described the observations that we will compare with the theoretical models, 
we can move on the expansion predictions for brane world gravity. 

3. Comparison of the DGP model vi^ith the data. 

In this section we will compare the predictions of the DGP model with the latest 
cosmological expansion history data set described in the previous section. As we have 
discussed in the introduction, the DGP model in its simplest form with one extra flat 
dimension ([ID gives rise to a modified Friedman equation as written in equation 

We have made a distinction between best model fit and parameter determination - 
when we talk about 'confidence levels' we are referring to the regions within which the 
values change from the minimum value by the critical amounts that are usually 
used (Xmm + 2-3 for 68% etc. ). If we use the term '1 cr', '2 cr' or '3 cr' we are referring to 
respectively the 68%, 95% and 99% significances obtained by comparing the to the 
number of degrees of freedom. If the errors were Gaussian and the model being fitted 
was a good one then those two measures of statistical significance should be equivalent. 
However, since we do not know if both those criteria are exactly fulfilled, it is important 
that the reader understands the distinction between the different quantities. 

3.1. Mathematical preliminaries 

It will be interesting to also consider generalisations of the DGP model which might 
result from having a higher dimensional bulk. Although such models have not been 
derived exphcitly, it is possible to guess at their possible form and the way that the 
Friedman equation would be modified as a function of the number of extra dimensions 
[23] . Actual realisations of higher dimensional models of this kind have problems which 




(5) 
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arise when one considers a delta function of stress energy at the position of the brane 
[22]. This is solved assuming a non-zero brane thickness, which means that a potential 
is required to maintain a nearly massive 4D graviton on the brane [61]. This potential 
will distort the spectra of gravitons in the extra dimensions and also therefore the 
leakage of gravity into the extra dimensions. It is not known precisely what form the 
modifications to the Friedman equations would take in such a situation in the cross over 
from 4D to higher dimensional physics, we assume that such a regularisation would take 
the generalised form [62] 

rr2 I ^ Pm /„x 

L\(3 + {LHY-^) 3Mj,i 
for zero spatial curvature on the brane. Here L corresponds to the crossover length 
scale and n is the number of extra dimensions. The extra term j3 parameterises the 
regularisation in cases where n 1. 
Dividing through by H^, we have 

^'(^) ^S-n) ... . = ^m(1 + Z? (7) 

where VLl = (HqL)"'^ and E{z) = H{z)/Hq. It is easy to see that when (5 dominates the 
denominatoii§|, the f2£,-term will behave very much like a cosmological constant term. 
As usual, we can constrain one of the parameters in terms of the others by setting z = 
in (I7j), giving 



P=^^-V^L (8) 

which is the equivalent of the equation f^A + = 1 in flat ACDM. 
3.2. Models with (3 = 0, including original DGP model 

The first models that we constrain are those with the regularisation parameter j3 set 
to zero. If there are too many parameters in a dark energy model it is usually rather 
easy to fit the data for some combination of those parameters, so one can assume priors 
like imposing a fiat universe in order to obtain interesting constraints. Here, because 
we have the same number of free parameter^ in these f3 = models as in the case of 
ACDM, we are able to include curvature and still get interesting constraints. We do this 
simply by replacing the Hubble constant squared with H"^ + fca~^ [13] . This gives us 




Pm 



(9) 



which for the case of the original DGP model leads to 

E\z) = ^k{1 + zf + [^^m{1 + zf + ^Il/A + VTvi)' (10) 
§ try saying that five times quickly 

II we have theoretical motivation for n being a discrete integer and each n corresponding to a different 
model 
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MLCS (162 data points) 


n=l 


ACDM (n=2) 


n=3 


n=4 


xLn(SNe) 


188 


188 


188 


188 


xLn(SNe+flat) 


188 


188 


189 


189 


xLn(SNe+flat+CMB) 


200 


188 


190 


197 


xLn(SNe+flat+CMB+BAO) 


215 


192 


191 


197 



Table 1. Best fit values for the fits of different data sets to tlie various models. SNe 
is the supernova data set (analysed using the MLCS2k2 procedure) without inclusion 
of the redshift error due to the peculiar velocities, next flatness is assumed then the 
CMB shift parameter is added to the data set. Finally we include the Baryon acoustic 
oscillation result (BAO) from the Sloan Digital Sky Survey. 

where the normal definition VL^ = ~k/{aHQ)'^ has been used. 

The different models were compared with the two sets of data, one of which 
contained supernovae that had been analysed using the SALT method of determining 
magnitudes while the other using the MLCS2k2 method. The data-set analysed by 
the ESSENCE collaboration with SALT gave rise to anomalously large values of 
seemingly due to a few outlying data points. Since the data analysed using MLCS gave 
per degree of freedom rather close to one we have more confidence in this data set 
and have restricted ourselves to using it. For this data set the intrinsic error in the 
magnitudes is taken to be 0.1 and the peculiar velocity error in the redshift, which 
has importance for the lowest redshift supernovae, is assumed to be 400 km/s. To 
demonstrate the effect of different error assumptions we present the resulting fits with 
and without inclusion of the peculiar velocity error. 

Figure [1] shows the comparison of the data (supernovae, galaxy survey baryon 
oscillation and CMB) with the DGP model and its variants where the supernova data 
has been analysed using the MLCS2k2 method. The confidence regions refiect the 
reported statistical uncertainties of the various measurements. The potential impact 
from systematic effects is not addressed by this analysis. 

Note that at this level of expansion history the n = 2 DGP variant is identical 
to ACDM although presumably perturbations would grow very differently in the two 
models. The actual values corresponding to the two different models DGP (n=l) and 
ACDM (n=2), as well as the higher dimensional generalisations for n=3 and n=4 are 
listed in the table [H and [2] . The effects of adding the velocity errors are of course lower 
minimum values as well as a small shift in the best fit parameter values, explaining 
the increased area of the supernova confidence regions and their shift in parameter space 
seen in figure [TJ 

All models can fit the supernova data alone equally well, and more or less equally 
can fit the data assuming fiatness. The requirement of having to fit the CMB data 
singles out ACDM as being slightly favoured over the other models, although the 
values are such that the DGP model is roughly 1-2 a disfavoured, depending on whether 
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DGP n=1 ACDy (n = 2) 




G,^ 0,6 0,8 1,0 (,:,2 OA 0,5 0,8 1,0 

0^ % 



Figure 1. Supernovae from ESSENCE, SNLS and nearby sample analysed with the 
MLCS2k2 method compared with models described by equation [9] with (3 = which 
include the basic DGP model (n = 1). Because ACDM is at this level identical 
to the n = 2 case we plot the different models in order of increasing n. The pink 
concentric ellipses correspond to the 68%, 95% and 99% confidence regions from the 
latest supernova data ([3]) when the peculiar velocity error is not yet included, the blue 
bands border the 99% confidence region for the baryon acoustic peak data ([4]) while 
the green region borders the 99% confidence region for the CMB shift parameter ([5]). 
The pink dotted line corresponds to the 99% confidence region for the supernova data 
when we have included the 400 km/s error in the redshift. The black line corresponds 
to spatially flat universes. 
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MLCS (162 data points) 


n=l 


ACDM (n=2) 
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xLn(SNe) 


160 


160 
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xLn(SNe+flat) 


161 


161 


161 
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xLn(SNe+flat+CMB) 


171 


161 


163 


170 


xLn(SNe+flat+CMB+BAO) 


180 


163 


164 


170 



Table 2. Same as tabled] except here the redshift error of 400 km/s has been included 
in the supernova data set. 
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Figure 2. Comparison between tiie results of fitting DGP and ACDM to tlie SNLS 
and ESSENCE supernova data set (filled in 68%, 95% and 99% confidence regions) 
and the Riess 07 Gold set (dotted lines). The solid black line corresponds to spatially 
flat universes. 



the peculiar velocity error has been included or not, which is not statistically significant 
enough to allow us to claim that the model is ruled out. 

If we add the baryon oscillation feature to the data we find that we do seem to 
be able to disfavour the DGP model at the 1.5-3 a level, again depending on what 
errors have been included, but as we have already stated in section [2l the reconstruction 
of the peak in the correlation function from redshift space depends upon assumptions 
which are not valid in cosmologies where the equation of state is varying over time. 
At some point in the future it may be interesting to investigate such systematic effects 
quantitatively. 

Note also that even if we from the last row of table [1] draw the conclusion that the 
DGP model is ruled out at a 3 a level, the ACDM model is then disfavoured at the 2 
a level. Looking instead at table [2] the values obtained for DGP and ACDM give 
goodness of fit of 16% and 46% respectively. 

Figure [2] shows a comparison between the results of fitting the DGP model and 
the ACDM model respectively to the Riess 07 gold set, used in [17], and the SNLS 
and ESSENCE data set, used in this paper. Clearly, there are inconsistencies between 
the two published sets of SN data and we have therefore not attempted to combine 
them. We see that for supernova data and a fiat prior only, the Riess data makes some 
distinction distinction between the two models while the SNLS data does not. It is on 
including the prior of the CMB shift parameter that the DGP and ACDM are equally 
well favoured by Riess data while the result in this paper is that DGP looks slightly less 
favoured than ACDM. 
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flat DGP /3 = flat DGP (3 = 




0.1 0.2 0.3 0.4 0.5 0.6 0.1 0.2 0.3 0.4 0.5 0.6 



Figure 3. Best fit values of n for /3 = and different values of the matter density 
JIm- Flatness is assumed and the priors of CMB and BAO included. On the left the 
supernova data used is analysed by MLCS2k2 and on the right the data analysed by 
SALT. 

3.3. Best fit values of n. 

An interesting exercise is to take the case /? = and allow the number of extra 
dimensions n to be a free, non-integer parameter. We then calculate which values 
of n fit the data best when we allow VLm to be a free parameter. 

Figure [3] shows us that the best fit region when /5 = occurs for values of n between 
1.5 and 3. The case n = 2 corresponds exactly to a cosmological constant (whether (3 = 
or not), so this tells us that without f3 the best fit to the data lies somewhere around 
the ACDM models, in agreement with the previous conclusions of [T3] although the 
preferred value of n may have increased slightly. 

3.4. Models with non-zero (3 

As stated earlier, for the higher dimensional generalisations of DGP we should include 
the parameter j3 which, while required for theoretical consistency, is an extra free 
parameter which renders the theory less predictive and makes it easier to fit the data. 
For this reason, for the higher dimensional cases, we will consider spatially flat universes 
with non zero {3. 

The original DGP model corresponds to n = 1 and in that case there is no need for 
a regularisation parameter (3. In terms of large scale geometry of space time the n = 2 
case is equivalent to ACDM with or without (3. We therefore only consider n = 3, 4, 5, 6. 

In all of these four models, the confidence region in the parameter space of VLm and 
VLl centers around VLm = 0.26. Taking this as our value for VLm-, we obtain a minimum 
value of VLl for each n by looking at the 95% confidence level away from the best fit, 
corresponding to a minimum value of j3 that increases for increasing n, as shown in table 

El 

In figure H] we also show how the different cosmological constraints cut into the 
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95% 


n = 3 


n = 4 


n = 5 


n = 6 
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(0.8) 


0.95 


1.15 


1.25 


/?> 





0.23 


0.74 


1.05 



Table 3. Minimum values of and /3 (via equation ^ with flM = 0.26) with 95 % 
confidence. For n=3, the minimum value of fi^ is restricted only through the cutoff 
where /3 goes negative. 

DGP(,S) n = 3, 6 




12 3 4 5 

Figure 4. Constraints on the higher dimensional versions of DGP taking into account 
the regularisation parameter (3. The filled in contours are the confidence levels for the 
combined for the n — 3 case (see note at beginning of section). The combined 
for the n = 6 case is also plotted for comparison using black lines showing that for this 
case, zero (3 is significantly disfavoured. 

parameter space of (3 and fi^ by plotting the confidence regions for two of the higher 
dimensional models. 

This shows that for a higher number of dimensions, the fit gets worse and one will 
have to increase the value of P in order to get closer to the cosmological constant case 
which fits the data well. As (3 becomes more dominant, the dark-energy term in eq. (I7l) 
approaches VLl/(3 and we can see from figure H] that the gradient approaches ~ 0.7 for 
larger /3, which agrees with the usual best-fit value of ^a. 

4. Summary and Conclusions 

Data from the Supernova legacy survey analysed with the SALT method has previously 
suggested that the DGP model is marginally disfavoured relative to ACDM [13]. At the 
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same time, a larger data set including recent data from the Hubble space telescope and 
analysed using the MLCS approach has led to the conclusion that the DGP model is in 
fact perfectly safe [17i| . 

In this work, we have looked at the SNLS data and new data from the ESSENCE 
collaboration analysed with the MLCS2k2 algorithm as reported in [21]. Combination of 
this data with the CMB constraint suggests that the DGP model is slightly disfavoured, 
and becomes more dis-favoured if one can treat the baryon acoustic peak as a valid data 
point. The magnitude of the change in the position of the acoustic peak in the galaxy 
correlation function when one moves from a background cosmology with a constant 
equation of state to a DGP universe depends both upon the different geometries and 
the way that structure grows in those universes. Since the growth of structure in a 
DGP universe may be rather different from in ACDM we are not able at this stage to 
say whether or not the DGP model seems to be marginally or significantly disfavoured 
using this data. Either way, we find that the tests of the DGP model yield significant 
differences when using the ESSENCE supernova data and the "gold set" [T71 116] . 

The fact that the two data sets lead to different conclusions about the same model 
is very interesting and outlines the challenges which need to be overcome in order to 
move into the era of precision Dark Energy measurements. 
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